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In the light-front formulation of field theory, it is possible 
to write down a chirally invariant mass term. It thus appears 
as if one could solve the species doubling problem on a light- 
front quantized transverse lattice in a chirally invariant way. 
However, upon introducing link fields and after renormaliz- 
ing, one finds exactly the same LF Hamiltonian as if one had 
started from the standard Wilson action in the first place. 
The (light-front) chirally invariant transverse lattice regular- 
ization is thus not chirally invariant in the conventional sense. 
As an application of the Wilson formulation for fermions on 
a _L lattice, we calculate spectrum, distribution functions and 
distribution amplitudes for mesons below 2GeV in a trun- 
cated Fock space. 



I. INTRODUCTION 

There are both phenomenological as well as theoretical 
reasons to study light-front (LF) quantization: degrees of 
freedom that are close to the relevant degrees of freedom 
in many high energy scattering experiments and the sim- 
plified vacuum structure. For recent reviews on these 
topics see Refs. (!]-§]• 

One approach that seems particularly promising to- 
wards performing numerical calculations of hadron 
bound states in QCD^ + i is the transverse (_L) lattice 



For our purposes, the most crucial difference between 
covariant field theories and LF field theories is the fact 
that LF field theories allow for [J] chirally invariant mass 
terms 
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In the canonical quantization approach, this term arises 
naturally in addition to the (chirally odd) helicity flip 
interaction term, which is also proportional to the mass. 

Since it is possible to write down a chirally invariant 
mass term, it is also possible to write down a chirally 
invariant analog of the Wilson r-term 



1 In fact, proper renormalization requires them Jl| 
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where oX nxi „ H is the lattice approximation to 
—ar'ipj^—dj_ip. One might thus hope that chiral fermions 
can be formulated on a transverse lattice without the 
notorious species doubling problem [||. The purpose of 
this paper is to give an explicit example which shows that 
this hope is premature. In fact, it will be shown that, 
after renormalization, one recovers the standard Wilson 
term in the interaction. 

In order to keep the complexity of the equations at a 
minimum, the lattice regularization is first explained on 
a _L lattice in 2+1 dimensions. This has the advantage 
that one does not need to specify the direction in which 
the transverse derivative acts. After extending the for- 
mulation to 3+1 dimensions, we proceed to apply Wilson 
fermions to a simple _L lattice model for mesons. 



II. FERMIONS ON A TRANSVERSE LATTICE 

A. Naive Fermions 

We start by rewriting the Lagrangian for a free fermion 
(2+1 dimensions, 4-spinors for the fermions) in terms of 

the LF-projections (f> = 7 2 7 ip and x = 7 2 7 V' 
C = i/j (i jd — m) vb 

= V2^id + <p + V^x^id-x - t (i(*xd x + (3m) x 
-xHia x d x + Pm)(j), (2.1) 

where a x = 7°7a; and (3 = 7 . Upon discretizing the 
transverse (x — > n) coordinate and replacing transverse 
derivatives by finite differences in Eq. (2.1) one thus 
finds the discretized action for (naive) free fermions on a 
transverse lattice 
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where (the dependence on the continuous longitudinal 
coordinates x ± is suppressed for notational convenience) 
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Since Eq. ( |2.3| ) does not involve a "time" {x + ) derivative 
of x, it is convenient to use the constraint equation 
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to eliminate \ prior to quantization, yielding the effective 
action for the dynamical (0) degrees of freedom 
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A simple plane wave ansatz shows that the above naive 
Lagrangian results in a dispersion relation which exhibits 
species doubling 
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where — — < k x < — . Thus, even for a — > 0, modes 
with k x f» 7r I a have finite energy. Of course, this has 
nothing to do with the LF approach chosen here, but 
is just the familiar phenomenon of species doubling for 
Dirac particles (first order derivatives!) on a lattice. 

The species doubling problem on a _L lattice has been 
dealt with in Ref. || using a generalization of the fa- 
miliar Kogut-Susskind fermions [ fio| |. However, in this 
paper, we will try different approaches, which are based 
on generalizations of Wilson fermions |ll[] . 



B. Ordinary Wilson Fermions 

The most familiar "patch" to deal with species dou- 
bling is to start from the naive action, plus a "Wilson 
r-term" of the form 

r - 

a 

= ~ {0n/3 [Xn+1 ~ 2X77 + Xn-l] 

a 

+X j n P [077+1 - 20„ + 0„_i]} 

where r = 0(1). The well known problem with the Wil- 
son r-term is that it has to have the Dirac structure of a 
mass term in order to eliminate "doublers" both for par- 
ticles and anti-particles and therefore chiral symmetry is 
explicitly broken. 

Including this Wilson r-term fl2.8| ), the constraint equa- 
tion for x reads 
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yielding for the effective Lagrangian of the dynamical 
degrees of freedom (0 n ) 
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The dispersion relation, with the r-term included, reads 
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and is, by construction, free of species doubling at k x 
±ir/a. 



C. Gauging Ordinary Wilson Fermions 

An action that reduces to the QCD coupling of 
fermions to transverse gluons in the naive continuum 
limit is obtained by introducing link fields U n , defined 
on the links that connect the site n with n+ 1, using the 
rule 



Tpn+llpn ► ^pn+lUllpn. 
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To keep the notation simple, we suppressed the depen- 
dence of U n on x~ . The coupling to the longitudinal 
gluon components is a separate issue and does not affect 
the main conclusions in this section. We will thus omit 
couplings to A ± here. Obviously, the resulting gauge in- 
variance includes only gauge transformations that do not 
depend on x ± . However, but if one wants to one can eas- 
ily remedy this by introducing couplings to A — which 
is what we will do in the model calculations in Section 

Hi- 
One thus finds for the action including link fields 
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After eliminating \n one thus finds 
where 

A//2 
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with A/ = m+ 2r/a, is the term containing no link fields 
and 
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is the term linear in the link fields. Once one renormal- 
izes, the term quadratic in the link fields [Cf^ 2 ] is deter- 
mined by the other two terms [pj[12"| and thus does not 
need to be specified for the point we are trying to make 
in this paper. 



D. A Modified r-Term 

In the context of LF quantization, it is natural not 
to write down the usual covariant Wilson r-term, but 
instead to modify only the kinetic mass term in the La- 
grangian after \ has been eliminated, i.e. instead of Eq. 
(2.8) we introduce a modified r-term, which is formulated 
directly in terms of the dynamical degrees of freedom 
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which is the discretized version of arib^ 
modified r-term yields a dispersion relation 
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which is very similar to Eq. (2.11), but nevertheless does 
not exhibit species doubling. What makes this "modified 
r-term" especially attractive is that it is chirally invari- 
ant (on the LF — see Appendix^ This feature can be 



best seen by re-expressing Eq. (2.18) in four component 
notation, yielding 
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At first, this seems to contradict the Nielsen- Ninomiya 
theorem jL|, which states (under very general assump- 
tions) that there is no chirally invariant way to solve the 
specie s dou bling problem. However, the additional term 
[Eq. ( |2.18 )] is non-local due to the inverse longitudinal 
derivative and thus the Nielsen-Ninomiya theorem does 
not apply. 



In order to render Eq. ( 2.18 ) gauge invariant, we re- 
place 
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and similarly for other non-local terms. Alt houg h the 
gauge invariant extension is not unique, Eq. (2.21) rep- 
resents the minimal extension of Eq. (2.18) that is gauge 
invariant ^, hermitian and invariant under transverse par- 
ity {n — > —n). Using this rule, we find the following terms 
contributing to the action based on this modified r-term 
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for the term containing no link fields as well as 
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for the term linear in the link field. As in the previous 
section, we do not exhibit the term quadratic in the link 
field since it follows uniquely in the renormalization pro- 
cedure. 

A comparison with the results from the previous sec- 
tion shows immediately that the naive r-term as well as 
the modified r-term yield, up to a redefinition of param- 
eters, identical expressions for C^'°, C^' 1 and C^'°, C^' 1 . 
Since these are just bare parameters, this means that 
up to and including linear terms in the link fields, the 



2 In fact, one can always start from an un-renormalized LF- 
Hamiltonian without "Compton terms" (i.e. the terms bilin- 
ear in both fermion and boson field) and then generate the 
Compton terms as counter terms. 



3 In the LF approach, after eliminating constrained degrees 
of freedom, already the canonical Lagrangian/Hamiltonian is 
non-local in the longitudinal direction. 
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independent gauge 
For the general case the same comments 
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standard and the modified Wilson terms yield the same 
interaction terms. Since the terms quadratic in the link 
fields are determined in the renormalization procedure, 
this means that the standard and modified Wilson for- 
mulations on a transverse lattice are identical. 

This is the central result of this paper, since it shows 
that once we introduce link fields and once we deter- 
mine the "Compton-terms" in the LF Hamiltonian in 
the renormalization procedure, also the terms quadratic 
in the link fields are the same for the naive Wilson r- 
term and the "modifie d" W ilson r-term. Even though 
our initial ansatz [Eq. ( 2 . 18))] looked very promising, all 



that has been achieved is that we reproduced the same 
interaction terms as the standard Wilson ansatz. Com- 
pared to the standard Wilson ansatz, nothing has been 
gained by trying to remove species doubling using the 
(LF-) chirally even kinetic mass term. 



step the anti-quark follows behind by absorbing the link 
field quantum. Propagation over more than one link pro- 
ceeds by repeating this process. The resulting sequence 
resembles the contracting and stretching motion of an 
inchworm. The similarity in the propagation mechanism 
and the fact that the typical length scale of the mesons 
is not one inch but one Fermi, motivated the name for 
the model. 

The kinetic part and self-interactions of the link fields 
have been introduced and discussed in Refs. ||]l4|] anc ^ 
will not be discussed here since the focus of this paper is 
on the quarks here. For details, the reader is encouraged 
to consult these References. 

As far as the longitudinal dynamics is concerned, in- 
troducing quarks leads to both a kinetic term 



quarks 
kin 



£4 



(3.1) 



E. From 2+1 to 3+1 Dimensions 

The main work in this paper so far was to show that 
the fermion terms up to and including linear in the link 
fields are the same for the standard and the modified 
Wilson formulation. Since terms nonlocal in both trans- 
verse directions {x and y) at the same time (x and y) 
involve 2 link fields, these do not occur when one focuses 
on terms up to and including linear. Therefore, the 3+1 
dimensional analysis turned out to closely parallel the 
2+1 dimensional analysis. Since no new physics insight 
was gained — in fact, the essential steps became more 
mysterious because of the more lengthy equations — the 
details of the 3+1 dimensional analysis were omitted in 
this work. The final result, i.e. equivalence between stan- 
dard Wilson and modified Wilson formulation, was the 
same as in 2+1 dimensions, since coupling constant co- 
herence between three point and four point couplings of 
bosons to fermions holds in both cases. 



III. THE FEMTOWORM MODEL FOR MESONS 

The coupling of quarks to the link-fields via the canon- 
ical terms leads to helicity flip (7-structure: a±), while 
the r-term does not flip the helicity of the quarks (7- 
structure: 0). In order to investigate the ramifications 
of these basic features for the interactions of quarks, let 
us consider mesons in the approximation that there is 
at most one link field separating a quark and an anti- 
quark. In this approximation, the quark and anti-quark 
are always on the same site or on neighboring sites and 
mesons propagate by one-boson exchange. For example, 
suppose that the quark and the anti-quark start out on 
the same site and that the quark hops to a neighboring 
site. Gauge invariance requires that the quark emits at 
the same time a link field quantum on the link that is 
now connecting the quark and anti-quark. In the next 



as well as a contribution from the quarks to the vector 
current 
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where i, 7 are color indices and where J J (U) is the 
vector current from the link field which has been in- 
troduced in Ref. §. The main difference between the 
longitudinal gauge field coupling of the quarks and the 
link-fields is that the quarks couple only to the longitu- 
dinal gauge field at one site, while the link fields couple 
to the longitudinal gauge field at the two sites that are 
adjacent to the link. In general, i.e. without truncation 
of the Fock space, the physical space of gauge invariant 
states thus consists of closed loops of link fields plus open 
strings of link fields with quarks and anti-quarks at the 
end. For studying mesons in the large Nc limit, we can 
focus our attention on quarks and anti-quarks connected 
by a string of link fields. The main approximation in this 
section will be to limit the length of this string to one 
lattice spacing. At first this seems to be a rather drastic 
approximation. However, if one considers that the lattice 
spacing in cont empor ary _L lattice calculations is on the 
order of 0.4/m Jujjljj], limiting the length of the string to 
one lattice unit seems to be a reasonable approximation 
for a first study of meson spectra on a _L lattice. 

Since the coupling of the quark fields to the longitudi- 
nal gauge field is local in the _L direction, the _L dynam- 
ics of the quarks is solely described by the hopping term. 
Using the result from Section || one finds the following 
for the Dirac part of the hopping terms: (Overall factors 
and the obligatory link-field operators are suppressed for 
notational convenience! It is implicitly understood that 
hopping is always accompanied by the emission or ab- 
sorption of a link field quantum on the link connecting 
the initial and final site of the hopping process) 
• naive term: 

hopping in x-direction: 
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-bl(ft±,Pf)bi(fi± + i x ,pf) 
-b\(n±,pj)bi(nx - i x ,pf) 
+b[(n±,pj)b^(n± - i x ,pf) 
-d\{n±,p'f)d l (ri ± + i x ,pf) 
+dl(n±,pj)d^(H± + i x ,pf) 
+d\(n ±) pj)di(nx - i x ,pf) 
-dl(n±,pj)di(nx - 4, Pi") 



hopping in y-direction: 
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+b\(ri_ L ,p+)b 1 (ri ± +T y ,pf) 
-b\(H±,pj)bi(n± ~i y ,pf) 
-b\(nx,pf)b^(fi± - i y ,pf) 
-d\(Hx,pj)di(n± +iy,pf) 
-d\(ri ± ,p+)d^(n± +iy,pf) 
+d](n±,p~f)d l (n ± - i y ,pf) 
+d t l (n ± ,p+)d- l (ri ± - i y ,pf) 



• r-term: 
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+6{(n_L,Pj)64.(nx + 4? Pi") 
+d\(n±,p^ )d 1 (n ± - i x ,pf) 
+d\(n±,p~f)d l (ri ± -i x ,pf) 
+b\{nx,Pf)bt(n± + T y ,pt) 
+b\(n±,p~f)b l (n_ L + i v ,pf) 
+d\(n_ L ,p+)d^(n ± - i y ,pf) 
+d\(n±,p~f)d l (n 1 _ - i y ,pf) 

Note that in the naive term, hopping in the the positive x 
or y direction (n — > n + 1) yields the opposite sign from 



hopping in the negative direction (n — ► n — 1). In the 
r-term, the sign is the same for hopping in all directions. 
This observation will become important when we discuss 
interference between various contributions. 



A. Solutions with P, 



= 0. 



Because of manifest (discrete) translational invariance 
in the _L direction, the total transverse momentum pj° ta/ 
is conserved. As we will see below, the case pt° tal = Q 
yields particularly simple bound state equations and in 
the following we will focus our attention on this case. 

The first simplification is a non-interference between 
naive hopping and r-term hopping: In general, one can- 
not distinguish between a link quantum emitted by the 
naive term and one that was emitted from the r-term. 
Thus a link quantum emitted from the naive term when 
the quark hops by one lattice unit may be absorbed 
by the r-term when the anti-quark follows behind and 
vice versa. Since the naive term flips the spin, but 
the r-term does not, this interference results in a non- 
conservation of the meson spin (which is of course a lat- 
tice artifact). However, this interference term has oppo- 
site signs for hopping into opposite directions. Therefore, 
for p*° ta/ — 0, there is an exact cancellation between in- 
terference terms when the meson hops in opposite direc- 
tions and there is not interference between link quanta 
emitted from the naive term and those emitted from the 
r-term. 

Note that, since the naive term flips quark helicity and 
the r-term does not, the above result has an important 
consequence for the spin structure of the qq (2 particle) 
Fock component: Because of the non-interference men- 
tioned above, the only conceivable mixing of 2-particle 
Fock components could be between states where q and 
q have both spins up or both spins down, i.e. "double 
helicity flip" . However, since the spin matrix clement 
for hopping in the x-direction with double helicity flip is 
(— — 1 (3.3), while th e sa me matrix element for 
y-hopping is (+1)( — 1) = —1 (3.4), both terms cancel for 
a plane wave with Pj_ = and there is no double helicity 
flip hopping. As a result, the total spin in the 2-particle 
Fock component is conserved. 

From a more general point of view, the above result 
about S z conservation in the 2-particle Fock component 
emerges since the transverse lattice is invariant under 
rotations by tt/2 around the z-axis, plane wave solu- 
tions with p t ° tal = must be eigenstates under rotations 
around the z-axis by n/2 (with eigenvalues ±1 and ±i, 
i.e. in a sense J z is conserved modulo 4). Since there is 
no orbital or gluon field angular momentum when the qq 
pair is on the same site this means that S z , which can 
assume only the values and ±1, must be conserved. 

Another important observation concerns the sign of the 
spin splittings between mesons. For P x = one finds that 
x-hopping is attractive in the pseudoscalar channel (spin 
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wave function: fj. — if), since for example 6(T — >J.) = —i 
and dQ— >|) = — i, i.e. together with another (— ) from 
the spin wave function and another sign from the energy 
denominator the total sign is (— ), i.e. attractive. The 
same result holds for y-hopping (kp = 0), since for ex- 
ample b(|— >I) = 1 and d(|— *|) = — 1. Likewise one can 
show that naive hopping yields a repulsive interaction for 
the vector channel (spin wave function: | J, + Jt), i.e. de- 
spite all approximations, the Weingarten inequality |l6| ] 
is satisfied for our model. 



B. Pure qllq States 

In the previous section we have seen that (discrete) ro- 
tational invariance yields a drastic simplification in the 
spin structure for p*° ta * = solutions to our model. 
Again for pt° tal — 0, we will demonstrate in this section 
that there are some particularly simple solutions which 
have a vanishing 2-particle Fock component. 

Mesons which mix with qq states at the same site must 
necessarily have total helicity ^ h — 0, ±1. However, even 
in the one link approximation, it is very easy to construct 
states which cannot mix with a qq state at the same site. 
The generic form of such states is 

\h = 2) = ^2 \b\(n± + ix)d\(n±) + £&j(nj_ + i v )d^(n±) 

- &j(n_L - i x )d\(n±) - ib\(H± - i y )d^(n±) |0) 

(3.6) 



and similarly for h = — 2 with the longitudinal wave 
function (fc + -dependence) factorized and with the link- 
field implicitly included. The angular dependence in Eq. 
(3.6) resembles that of a state with L z = +1, but strictly 
speaking on a lattice there is of course no rotational sym- 
metry and hence no angular momentum . Ne vertheless, a 
state with orbital structure as in Eq. (3.6) cannot mix 



with states where a qq pair is at the same transverse site: 
first we note that the matrix element for hopping due to 
the r-term is independent of the direction and thus there 
is destructive interference when adding up the phases 
from the state for the quark to hop from different direc- 
tions to the site of the anti-quark (1 + £ — 1 — i= 0). The 
phase of the matrix elements of the naive term for a quark 
with positive helicity to hop to the site of the anti-quark 
are —i, 1, i, —1 for hopping from the +x, +y, —x, — y di- 
rection respectively. Together with the phase of the wave 



function, one thus obtains destructive interference for the 
amplitude due to the r-term as well 

(qq\H q _ hop \h = 2) cx HX+l) + (+l)(+£) 

+(+£)(-!) + (-!)(-£) =0. (3.7) 

Similarly, one can verify that the total amplit ude f or anti- 
quark hopping is zero for the h = 2 states (3J3). Note 
that it was crucial for the cancellation argument that the 
helicity of both q and q is positive. 

There are more states that do not mix with qq states 
for Pj° tal = 0. For example, any state with L z = 2, 
regardless of spin, i.e. any state of the form 

\L Z = 2)=J2 [ b Hn± + 4)rf f (nx) - b^(n ± + i y )d^(n ± ) 

+tf(n± - i x )d^(n±) - 6+(n± - i y )S(n^ |0) 

(3.8) 

has a vanishing matrix element with states where the qq 
state is on the same site. To see this, let us first consider 
r-term hopping. Since there are no phases for hopping 
in different directions from this term in the Hamiltonian, 
the only phases arise from the state and the add up to 
zero 1 — 1 + 1 — 1 = 0. This argument applies regardless 
of spin. The matr ix elements of the naive hopping term 
between Eq. (3.8) and a pure qq state vanish, since the 
naive hopping term in the Hamiltonian has always op- 
posite signs for hopping in opposite directions, the wave- 
function yields the same sign for separations in opposite 
directions. The total matrix element of the naive hop- 
ping term is thus zero due to cancellation from opposite 
directions. The pure qGq states are summarized in Table 
I. 

TABLE I. Quantum numbers of states that cannot mix 
with qq states. The C and P quantum numbers contain all 
signs from intrinsic parities, spin and transverse symmetry 
but do not include signs from the longitudinal orbital wave 
function. 



helicity 


L z 


S 


s z 


c 


p 


2 


1 


1 


1 


+ 


+ 


2 


2 


1 









2 


2 








+ 




1 


2 


1 


-1 







5 We define total helicity through the eigenvalues of the dis- 
crete rotation operator for rotations by n/2 around the z-axis. 
Since the resulting definition is unique only up to four times 
an integer, we break the ambiguity by always selecting the 
solution with \h\ < 2. 
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The total helicity of these states is obtained from the 
sum of the discrete orbital angular momentum and the 
spin (which is conserved modulo 4). The parity and C- 
parity quantum numbers given in Table I include all in- 
trinsic parities as well as signs from the spin and trans- 
verse wave functions. For a state with no nodes in the 
longitudinal wave function, the P and C parities given 
in the table are thus immediately equal to the P and 
C parity of the corresponding state. However, for more 
complicated longitudinal wave functions, the P and C 
parities in the table need to be multiplied by the signs 
from the longitudinal wave function. 

The longitudinal dynamics of all the states in Table I 
are described by the same relatively simple wave equa- 
tion, which reflects the fact that (note Nc — > oo) the 
quark and anti-quark interact only with the link field via 
the 1+1 dimensional gauge interaction 



M 2 ^(x,y) 



ml - G 2 ml - G 2 



1 — x — y 



ip(x,y) 



G 



1 v dz (1 — x — y) + (1 — z — y) 



(x - z) 2 2^/l - x - y^/1 - z-y 
2 f 1 -* dz (1 - x - y) + (1 - z - y) 

where x, y are the longitudinal momentum fractions of 
the quark and anti-quark respectively and where G 2 — 
The integrals in Eq. (3.9) are principal value inte- 
grals. 

Note that the link fields should not be directly inter- 
preted as gluons and therefore the mesons in Table I, i.e. 
the solutions to Eq. (3.9), are not necessarily hybrid or 



exotic mesons even though they cannot mix with pure qq 
states. The physical reason why these states cannot mix 
with qq states on the transverse lattice is that conserva- 
tion laws (e.g. angular momentum) do not allow that q 
and q sit on the same transverse coordinate, and gauge 
invariance on the _L lattice requires the presence of link 
fields. However , hy brid mesons correspond to those so- 
lutions to Eq. (3.9) where the link field is in an excited 
state. 



C. Numerical Calculations 

In the numerical calculations, we solved the cou- 
pled Fock space equations, including qq and qUq states 
projected onto vanishing total transverse momentum 

ptotal q 

Furthermore, after making this plane wave ansatz, one 
can reduce the equations of motion to equations that 
resemble the coupled Fock space equations in collinear 
QCD (with truncation to 3 particles) except that the 
link fields not only have vector couplings to the quarks 
but also scalar couplings from the r-term. Since collinear 



QCD models have been studied extensively in the litera- 
ture, we omit these details here and refer the interested 
reader to Refs. @@. 

In the Hamiltonian, we dropped the seagull interac- 
tions involving instantaneous fermion exchange! As has 
been explained in Ref. Jl9|j20|l , gluon loops connected to 
gluon emission and absorption vertices are important for 
dynamical vertex mass generation. If one works with 
a truncated Fock space, such vertex corrections are in- 
cluded in an asymmetric way: for example, if one limits 
the Fock space to qq and qll q then vertex corrections to 
boson emission can only be generated by the Hamilto- 
nian in the "in" state but not in the "out" state because 
the "out" state contains already the emitted boson (and 
conversely for absorption). As a result, in and out states 
are treated asymmetrically. In order to understand the 
consequences of such an asymmetric treatment, consider 
the momentum dependence of the bare emission vertex 
(helicity flip) 



T. 



emission 



X 



iiq 




(3.10) 



Loop corrections in the initial state connected to the ver- 
tex by instantaneous interactions renormalize only the 
mass term which gets divided by pf but not the mass 
term divided by p~f (and vice versa for loop corrections 
in the final state) |2C|] . However, this is exactly what 
happens to an emission vertex when the Hamiltonian 
is diagonalized in a Fock space truncated basis. There 
are several possible remedies one might think of, such as 
perturbatively including loop corrections which are sup- 
pressed by the Fock space truncation. |^| However, in this 
paper we chose the most simple solution and dropped 
instantaneous fermion interactions altogether, which re- 
moves loop corrections connected to vertices by instan- 
taneous interactions completely. Since we are aware that 
those interactions are crucial for dynamical vertex mass 
generation [ ft9|]20| , we compensate for this omission by 
using constituent quark masses (300MeV) for the vertex 
mass and not current masses. Note that since there is 
also an unspecified vertex coupling constant multiplying 
the link filed emission terms, the precise numerical value 
of the vertex mass is not crucial. However, it is impor- 
tant to note that we used a value that does not vanish 
in the chiral limit. Furthermore, it is important that, for 
simplicity, we do not allow this coupling to run |pof . 

Furthermore, in the qUq Fock component, we use a 
constituent mass for the kinetic mass of the quarks as 
well. In the qq component of the Fock space, where one 
has to add an infinite counter-term for the one loop self 



6 Note that DLCQ does not completely solve this problem, 
even if all Fock components allowed by the discretization are 
included, because of the loss of boost invariance in subgraphs. 
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energy, we use the finite part of the kinetic mass as a fit 
parameter to fit the 7r and p masses. 

In general, the self interactions of the link field are 
described by a complicated effective potential [O . How- 
ever, in the one link approximation only the quadratic 
term (the mass term) contributes, which we take to be 
about 300MeV as well, which lies in the range of val- 
ues suggested in Ref. Q. In the spirit of an effective 
link field interaction (obtained in principle by "smearing" 
and integrating out short range fluctuations) it would be 
consistent to also introduce effective interactions for the 
quark fields. For example, an effective interaction in- 
duced by short range fluctuations could be a spin depen- 
dent "contact interaction" similar to Gross-Neveu inter- 
actions. It is also conceivable that one could introduce a 
spin dependent contact term such that it partially can- 
cels the violations of rotational invariance. However, in 
this first study case, we wanted to explore how well one 
can reproduce the light meson spectrum with canonical 
terms only and we considered canonical quark terms only. 

The most important scale in _L lattice calculations 
is set by the dimensionful coupling of the quark and 
link fields to the instantaneous Coulomb interaction act- 
ing on each site. If one suppresses higher Fock com- 
ponents then this coupling directly translates into the 
longitudinal string tension (a = ^G 2 , which we fix at 
(7 = IGeV/ fm. We have also tried smaller values of the 
string tension, but then we were forced to use a very 
large coupling for the hopping term in order to generate 
the experimentally observed ir — p splitting and the spec- 
trum became numerically unstable. This is because, for 
a given quark mass, smaller string tensions lead to wave 
functions which vanish faster at the kinematic end-points 
and hence lead to smaller contributions from the one bo- 
son exchange interaction which is also singular at these 
end-points. 

In our numerical calculations we applied a sharp cutoff 
to the matrix element for hopping (i.e. link field emission 
and absorption) processes if the momentum of the final 
state quark is smaller than a fraction e of the momentum 
of the initial state quark. We varied this cutoff between 
e = 0.2 and e = 0.1 and verified that the spectra are sta- 
ble with respect to variations of this momentum fraction 
cutoff. 

As long as the r-term has a coupling which is about 1 /4 
of the coupling for the naive hopping term or smaller, the 
low energy spectrum for P± = is rather insensitive to 
its precise value. In fact, for such values, the main effect 
of the r-term is to raise the energy of mesons consisting 
of "doubler quarks" from the edge of the Brillouin zone. 
When the r-term was taken to be larger (about 1/2 rel- 
ative to the naive hopping term or more) there was also 
a strong effect on the low energy spectrum for P± = 
and we were no longer able to simultaneously fit the ir 
and p meson masses. We thus performed the numerical 
calculations where the value of the r-term coefficient is 
1/4 of the naive term coefficient. 



Given all constraints listed above, we fitted the only 
remaining free parameters — the finite part of the ki- 
netic mass counter-term and the vertex mass times the 
coupling of quarks to A± — to the ir mass as well as the 
p mass. Since rotational invariance and parity are not 
manifest symmetries, it is necessary to be more specific 
here about what we mean by the tt and the p mesons. The 
7r-meson is taken to be the lightest meson with C = + 
and helicity h — 0. Since the p meson has spin 1, it can 
have both helicity h — as well as h = ±1. Since our ap- 
proach as well as model approximations break rotational 
invariance, we do not find these states to be degenerate: 
the lightest meson with h = 1 and C = — and the lightest 
meson with h = and C = — are not exactly degenerate. 
In the fitting procedure, we chose to fit the center of mass 
of these two solutions to the experimentally observed p 
meson mass m p = 770MeV. 

The resulting meson spectrum is shown in Fig. [l] for 
mesons with positive C-parity and in Fig. || for mesons 
with negative C-parity. 
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FIG. 1. Comparison of the calculated meson spectrum for 
helicities h = 0, 1, 2 with the experimental meson spectrum for 
states with positive C-parity. Those states where an obvious 
assignment for the parity quantum number seemed possible 
are labelled with a '+' or while states with ambiguous 
parity quantum number are labels with a '?'. Degenerate 
states are displaced slightly in the vertical direction in order 
to make the degeneracy visible. 

Unlike C-parity, the usual parity is not a manifest sym- 
metry in the LF framework. For states with a large 
2-particle component, we used the symmetry of the 2- 
particle component to assign a parity quantum number: 
for example, let us consider the vector current, which can 
obviously only couple to negative parity states. The vac- 
uum to meson matrix elements of the good component 
-07 + -0 of the vector current operator is proportional to 
J dx(f) n {x), where 4> n {x) is the 2-particle component of 
the wave function. Obviously, this matrix element can 
be non- vanishing only if the 2-particle component of the 
wave function is even under x <-> (1 — x). Similar obser- 
vations can be made for positive C-parity operators. We 
therefore assign a negative parity to states with an even 
wave function in the 2-particle component of the Fock 
space, and a positive parity to states with an odd wave 
function in the 2-particle component. Note that this as- 
signment is also consistent with simple quark model con- 
siderations (i.e. free particle parity). 
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FIG. 2. Same as Fig. |y, but for negative C-parity. 

For states with very small or vanishing 2-particle com- 
ponent, this procedure becomes dubious or useless. In 
those cases we determine the parity of a state from the 
symmetry of the 3-particle wave function: First there is 
transverse parity (which is a manifest symmetry both in 
LF quantization and on the _L lattice). Secondly, for the 
"longitudinal parity" we combine the intrinsic parity of 
a qq pair (— ) |^] with the orbital parity. The orbital par- 
ity is determined using the free (no interactions) parity 
operator. Note that instead of applying the free parity 
operator, one can simply count nodal lines: in addition 
to transverse and intrinsic parity, there is a (— ) if the 
3-particle wave function is odd under qq exchange and 
there is another (— ) if there is an odd number of nodes 
in the wavefunction describing the motion of the link field 
relative to the "center of mass" of the qq pair. For ex- 
ample, a state with positive transverse parity and where 
the longitudinal qU q wave function has no node at all 
would be assigned P = — . Of course, counting of nodes 



7 Note that there is no intrinsic parity for the link field. 
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is sometimes ambiguous, but in those cases using the free 
parity operator is also not very conclusive. 

However, while the above parity assignments proce- 
dure seemed reasonable to us, it is not completely un- 
ambiguous and we thus plotted states with positive and 
negative parities in the same figure. 

There are several observations one can make from Figs. 
|l| and |2[ First of all, because of the violations of ro- 
tational invariance, states that one would expect to be 
degenerate (such as the h=0 and h=l states of the p 
meson) are not degenerate. However, considering the 
crudeness of the model, the violations of rotational in- 
variance are actually moderate. The second observation 
one can make is that even though our model properly 
describes the gross features of the spectrum, there are 
clearly states "missing" in our model and the energies of 
excited states typically lie too high. The missing states 
can be easily understood by noting that the truncation 
of the Fock space to 3 particles suppresses degrees of 
freedom associated with higher Fock components. The 
result that excited states lie systematically too high also 
has to do with the fact that a lot of physics has been left 
out by our truncation of the Fock space. For example, 
restricting the maximum separation between quark and 
anti-quark to one link, has a stronger effect on excited 
states since those tend to be larger (if one allows them 
to spread out). 

Encouraged by the numerical spectrum, we proceeded 
to calculate other hadronic observables — particularly 
those where the use of LF field theory is advantageous: 
LF wave functions and parton distribution functions. LF 
wave functions or distribution functions for the 7r and 
the two helicity states for the p are shown in Fig. pi 



8 These are defined to be the longitudinal momentum space 
LF wave function (where momenta are measured in units of 
the total momentum) in the qq Fock component with zero 
transverse separation, i.e in our calculation the valence wave 
function. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.0 

p(h = 0) x c 




0.4 0.5 0.6 

P(h=1) x q 

FIG. 3. LF wave functions <j>2 (x q ) in the 2 particle compo- 
nent of the Fock space versus the momentum fraction carried 
by the quark x q for the tv, and the h = 0, 1 states of the p 
meson. 

What one observes in Fig. || is that all of them 
are rather flat with a maximum in the middle, which 
reflects the strong binding of these states. Thus de- 
spite the low scale (recall our large transverse lattice 
spacing a w .4/m) the pion wave function (^^(x) has 
a shape which already resembles the asymptotic form 
<pn{x) oc x(l — x), but it is much more flat. A closer 
look at the wave functions shows that the p meson wave 
functions are slightly more peaked in the middle than 4>^ 
which reflects the weaker binding of the p meson com- 
pared to the 7T. 

In Figs. |] and || we show the longitudinal wave func- 
tions of the 7T and the p in the 3-particle Fock compo- 
nents. Fig. ^ shows the component of the 3-particle Fock 
component where the quark or the anti-quark helicity has 
flipped compared to the valence component (generated 
by naive hopping!), while Fig. (|^) shows the component 
where the qq helicity wave function is the same as in the 
valence state. 
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around the z-axis (see the discussion in Section [II B) but 
with the glue in its ground state. The mass of this state 
is about 1.63 GeV, i.e. much higher than experimentally 
observed mesons with J = 2, which we attribute to the 
one link approximation. The second excited solution to 
the qU q equation (lower left in Fig. ||) has a node in the 
wave function for the link field relative to the qq pair, i.e. 
this state should be interpreted as a hybrid excitation. Its 
mass is about 2.27GeV, but again we expect significant 
corrections from higher Fock components which we have 
not included. Finally, the first and third excited state 
(upper and lower right in Fig. ^ respectively) correspond 
to states where there are nodes in the qq relative wave 
function, but not relative to the link field, i.e. we inter- 
pret these solutions again as non- hybrid excited states. 
Their masses are 2.01GeV and 2.32GeV respectively. 



p(h=1,t40 



p(h = 1 ,U) 



FIG. 4. LF wave functions for the 7r, and the h = 0, 1 states 
of the p meson in the Fock component with qq and one link 
field U and with helicity for the quarks flipped compared to 
the valence component. The x-axis is the momentum fraction 
xv carried by the link-field and the y-axis is the difference 
between the momentum fraction carried by the quark and 
the anti-quark x q — Xq. For the p meson with h = 1 there are 
two possible spin configurations for the quarks, which have 
different wave functions. 
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Ground State 






Excited State 



0.2 0.4 0.6 0.8 

3. Excited State 



p(h = 0) 



P (h = 



FIG. 5. Same as Fig. but for the 3-particle component 
with the same quark spin wave function as the valence com- 
ponent. 



In Section |IIIB| wc demonstrated that there arc states 
in the spectrum which do not mix with the qq component 
of the Fock space. In Fig. || we show the (longitudinal) 
LF w ave functions for the four lightes t so lutions to Eq. 
(3.9). The lightest solution to Eq. (3.9), displayed in 
the upper left of Fig. ||, has no nodes, i.e. we interpret 
this solution as a qq state with orbital angular momentum 



FIG. 6. LF wave functions for states that have vanishing 
qq Fock component. 

Quark distribution functions ("structure functions") 
for the 7r and p mesons are shown in Fig. 0. For the p with 
helicity h = 1, the polarized quark distribution function 
is also displayed. Note that, since the p meson has spin 1, 
even its unpolarized quark distribution function does not 
need to be the same for the h=0 and h=l states, i.e. the 
fact that our numerical results for the unpolarized dis- 
tribution functions of the p states with h = and h = 1 
should not necessarily be interpreted as a manifestation 
of rotational invariance violations! We should emphasize 
that the quark distribution functions should not vanish 
at the origin, because the interaction coupling the 2 par- 
ticle and 3 particle Fock components is singular enough 
as the fermion momentum goes to zero to overcome the 
suppression due to the rise of the kinetic energy in the 
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3 particle component when the momentum of one of the 
particles vanishes. However, in the numerical results this 
singularity is suppressed due to our cutoff on the mo- 
mentum fractions in a hopping (i.e. link field emission) 
process. Varying this cutoff e hardly changes the struc- 
ture functions except in the extreme vicinity of m = 0, 
where the "gap" fills up as e is lowered. The structure 
functions displayed are for a value of e = 0.1. The model 
yields for the momentum fraction carried by the quarks 
plus anti-quarks in the tt about 75.4 %, while the num- 
bers for the p are slightly higher with (x q ) p ( h=a ) » 77.4% 
and (x q )ptf l =i) ~ 76.2%, which again reflects the stronger 
binding in the n (resulting in a "smaller" object, where 
less volume is filled with gluon fields). In deep inelastic 
scattering experiments, the structure function of the it 
is usually determined at rather large values of Q 2 . if we 
were to evolve our distribution functions to larger values 
of Q 2 we would obtain momentum fractions carried by 
the quarks which are comparable to the experimentally 
measured value of about 50%. However, since the per- 
turbative limit of the _L lattice is still poorly understood 
E3 , we will not attempt to evolve our numerical results 
to higher values of Q 2 . 




0.2 0.4 0.5 0.S 1 0.2 0.4 0.5 0.8 

n x ^ p(h = 0) 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.S 



p(h= 1 ) x '< p(h = 1) 

FIG. 7. Unpolarized quark distribution functions for the 7r, 
the p with helicity h = and the p with h = 1. For the p 
with h = 1 the polarized quark distribution is also shown. 

In the polarized quark distribution function for the p, 
one can notice a sign change for small x. The reason for 
this sign change is the following. At large x the valence 
quarks dominate, which have the same helicity as the par- 
ent p). At small x, the distribution function is dominated 
by quarks (and anti-quarks) which have just emitted a 
link field and link filed emission via the naive term leads 



to helicity flip for the quark (or anti-quark). Since the 
coupling of the naive hopping term is taken to be larger 
than the coupling of the r-term, helicity flip dominates 
at small x over non-flip and thus the polarized parton 
distribution changes sign. In order to interpret this re- 
sult, we emphasize again that the parton distributions 
calculated should be interpreted at a very low momen- 
tum scale. Perturbative evolution would clearly cover 
up this sign reversal. However, it is interesting to note 
that, based on the above results, the spin fraction carried 
by the quarks spin in the p meson turns out to be about 
52.9% . Because of the large Nq approximation employed 
in this model, this number should be interpreted as the 
non-singlet spin fraction carried by the quarks. Unfortu- 
nately, no experimental data on the structure function of 
the p meson is available, but a number which has a sim- 
ilar physical interpretation would be the ratio between 
gA in the nucleon and its naive quark model value, which 
turns out to be « 0.77. We believe that the rela- 
tively small fraction of spin carried by the quark spins is 
an artifact of the constant (i.e. non-running) dynamical 
vertex mass which we have used for simplicity and we ex- 
pect that a more refined treatment leads to less negative 
polarization of the quarks at small x. 

IV. SUMMARY 

One of the main message of this paper is that even 
though it is possible to write down a mass term that is 
invariant under chiral transformations (on the LF), this 
does not help to solve the species doubling problem in a 
way that is chirally invariant in the usual sense. |^| 

We used the kinetic mass term, which is chirally in- 
variant in the LF approach, to remove species doubling. 
However the relevant 2 boson - 2 fermion coupling con- 
stants are not independent, but rather must be chosen 
so that they cancel divergences arising from iterating the 
three point interactions. Once one uses this result, the 
LF Hamiltonian based on standard Wilson fermions and 
the LF Hamiltonian based on a Wilsonian prescription 
(applied to the kinetic mass only) become identical. 

As a first application for the transverse lattice formula- 
tion of Wilson fermions, we studied the meson spectrum 
in a simple model where we truncated the Fock space to 
qq states (q and q on the same site) and qU q states (q and 
q on adjacent sites) In this approximation, we are able 
to separate the transverse motion, and, for pj° tal = 0, 
we derive an effective collinear equation for the model 
which can be easily solved numerically. After fitting the 
7r and p masses as well as the longitudinal string tension, 



9 Once again, since LF Hamiltonians are nonlocal in the lon- 
gitudinal direction, the Nielsen-Ninomiya theorem [^j does 
not apply directly. 
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we obtain a qualitatively reasonable meson spectrum, al- 
though excitation energies tend to be systematically too 
high. For the lowest hybrid meson state we obtain a mass 
of 2.27GeV, which we believe is too high due to our Fock 
space truncation. The distribution amplitudes for the it 
and p mesons are both peaked at x = 0.5 and vanishing 
at the end-points, with the p meson distribution function 
being slightly more peaked. The momentum fraction car- 
ried by the quarks (and anti-quarks) in the 7r meson is 
about 75% with its value for the p meson being slightly 
higher. For the p meson with helicity h — 1 we find that 
the spin fraction due to the spins of the quarks is about 
53% . 

Considering the simplicity of the model calculations, 
the numerical results on hadron observables are very 
promising and encouraging. However, because of the sim- 
plistic nature of our first study calculations, there are 
many improvements possible. First of all, because of the 
large lattice spacing employed in current _L lattice cal- 
culations, one might also consider introducing "improve- 
ment" terms jTJJ for quark fields. Secondly, especially 
for studying excited mesons, it would be very useful to 
include states with more than one link field in the Fock 
expansion, thus allowing the mesons to "expand" in the 
transverse direction. In connection with including higher 
Fock components, one also needs to address the question 
of dynamical vertex mass generation pp|] and the run- 
ning of the effective vertex mass in more detail in order 
to properly incorporate the physics of dynamical chiral 
symmetry breaking into the effective LF Hamiltonian for 
QCD on a _L lattice. 



where iD± = id± + gA± . Once tp- has been eliminated, 
chiral transforma tion s are applied to ip+ only! The con- 
straint equation (Al) is in general inconsistent with the 
usual chiral transformation, since 



75 



1 

Jd+ 



a±iD± + 7 to ip+ —— - d±iD± + 7 to 75^+ 
J iO + L J 

(A2) 

This explains how it is possible that the standard Wilson 
r-term, which breaks the usual chiral symmetry, does not 
break LF chiral symmetry. 

It is interesting to note that only for to = do chiral 
transformations on the LF and those defined covariantly 
agree with one another. In fact, it turns out that the 
generator for chiral transformations on the LF is noth- 
ing but the helicity of the fcrmions. This makes sense 
since for m — helicity and chirality become the same 
(as do chiral transformations defined covariantly and on 
the LF) and both become conserved. Given a LF Hamil- 
tonian, this fact makes it very easy to verify whether it is 
chirally invariant: if quark helicity is conserved by each 
interaction then the Hamiltonian is chirally invariant. 

Manifest (LF-) chiral symmetry implies that the helic- 
ity of each quark is conserved. This fact gives rise to pe- 
culiar degeneracies: for example, mesons with S z — and 
C = ±1 (the 7r and the S z = polarization states of the 
p) are exactly degenerate — a result which is phenomeno- 
logically not very desirable. Simple estimates show it is 
conceivable that a finite n— p splitting remains in the chi- 
ral limit due to the singular end-point behavior of wave 
functions and matrix elements of the relevant interac- 
tions. However, this point requires further clarification. 
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APPENDIX A: CHIRAL SYMMETRY ON THE 

LF 

Chiral symmetry on the LF and chiral symmetry in 
the covariant approach are not exactly the same (here we 
closely follow Ref. pi). In LF field theory, it is convenient 
to decompose ip = ip+ +"0-, where ip± = |7°7 ± V' an d 
to eliminate ip- using a constraint equation 



1>- 



1 



a±_ ( id± + gA± ) + 7°to 



+ • 



(Al) 



10 In the main text, we used 
the use of subscripts. 



4> = ip+ and x = 0- to avoid 
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